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ABSTRACT 

Aims. We study the influence of the environment on the evolution of galaxies by investigating the luminosity function (LF) of galaxies 
of different morphological types and colours at different environmental density levels. 

Methods. We construct the LFs separately for galaxies of different morphology (spiral and elliptical) and of different colours (red 
and blue) using data from the Sloan Digital Sky Survey (SDSS), correcting the luminosities for the intrinsic absorption. We use the 
global luminosity density field to define different environments, and analyse the environmental dependence of galaxy morphology 
and colour The smoothed bootstrap method is used to calculate confidence regions of the derived luminosity functions. 
Results. We find a strong environmental dependency for the LF of elliptical galaxies. The LF of spiral galaxies is almost environment 
independent, suggesting that spiral galaxy formation mechanisms are similar in different environments. Absorption by the intrinsic 
dust influences the bright-end of the LF of spiral galaxies. After attenuation correction, the brightest spiral galaxies are still about 
0.5 mag less luminous than the brightest elliptical galaxies, except in the least dense environment, where spiral galaxies dominate the 
LF at every luminosity. Despite the extent of the SDSS survey, the influence of single rich superclusters is present in the galactic LF 
of the densest environment. 



Key words, cosmology: observations ^ 
- dust, extinction 



1. Introduction 



cosmology: large-scale structure of universe ~ galaxies: luminosity function, mass function 



Understanding the formation and evolution of galaxies is one 
of the biggest challenges of observational cosmology. The lu- 
, minosity function (LF) is in this respect one of the most funda- 
mental of all cosmological observables, helping us to describe 
the global properties of galaxy populations and to study the evo- 
. lution of galaxies. The dependence of the LF on cosmic time, 
' galaxy type, and environmental properties gives insight into the 
physical processes that govern the assembly of the stellar content 
of galaxies. 

The first determinations o f the galaxy LF were made sev- 
. eral d ecades ago ( Kiang.,196l1;IChristensen|[T975l;lKirshner et all 
' Il979l) : in the following studies, the number of g alaxies used 
' to calculate the LF has increased continuously dTullvl Il988t 
Efstathiou et al. 1988; Loveday et al. 1992). The Las Campanas 
Redshift Survey measur ed the general LF of galaxies with a 
, good accuracy (iLin et al][T996l: [Bromlev etal.lll998t IChrisdeinI 
l2000h . 

Our current understan ding of the general LF owes much 



to the 2dFG RS (Norber g etal 
jBlanton et al.l |2003b 



20021) and SDSS 



surveys 

M ontero-Dorta & Pradal l2009l) . These 



new samples of galaxies make it possible to study the depen- 
dence of the LF on a large number of different galaxy proper- 
ties, as galaxy morphology, colours, star formation rate, local 
and global density environment, etc. 

The morphology of a given galaxy is a reflection of its 
merger history. In studies of the LF, galaxy morphology has 
been determined either by its co lours ( Yana et al. 2009), spectra 
dFolkes et al.l 1 1 9991: iMadgwick e t al. 2002; de Lapparent et al. 
12003 ) or photometric profile (iBell et al.i i2003i: iDriver et al.l 



| 2007ab ; some studies us e artificial neural networks for mor- 
phological classification (iBall et al.l 12006). The most accurate, 
but by far the most time consuming approach is to use visual 
classification (Marzke et al. 1994, 1998; Kochanek et al. 20Qll 
Cuesta-Bolao & Serna 2003; Nakamura et al. 2003). For the 
SDSS survey, visual classifi cation has become possible thanks 
to flie Galaxy Zoo project (iLintott et al.l l2008h that will help 
us to study the morphology and t he LF in detail in the future. 
Recently Huertas-Companv et al.l (1201 ih have published an au- 
tomated morphological classification based on the Galaxy Zoo 
data. In all these studies, the classification of early-type and late- 
type galaxies is different, but ah studies agree that later-type 
galaxies have a fainter characteristic magnitude and a steeper 
faint-end slope of the LF. The biggest differences in previous 
studies are found at the faint-end of the LF, where classification 
is less certain than for brighter galaxies. 

To understand how galaxies form, we also need to under- 
stand where galaxies are located; it is essential to study the 
LF dependence on the environment. It is well known from 
the halo occupation distribution models that t he local environ- 
ment is crucial for the galaxy distribution (e.g. Kandivarez et al.l 
2006; Park et al.. 2007): luminous galaxies tend to occupy high 
mass haloes and low luminosity galaxies reside mainly in low 
mass haloes. This motivates the study of the L F in galaxy 
groups (Xia et al . 2006; Zandivarez et al. 2006; H ansen et al.l 
120091: lYang et'ani2009h . 

A likewise important, but not so well understood factor is 
the global environment where the gal axy is located - its p lace 
in the supercluster-void network. In iTempel et al.l (l2009b we 
have found that the global environment has an important role 
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in determining galaxy properties. Some studies have bee n dedi- 
cated only to special regions: e.g.. lMercurio et al ] (12006') inves- 
tigate the Shapley supercluster. The dep endence o n the environ- 
ment has been also studied n umerical ly (Mo et al. 2004) and us- 
ing se mi-analytical models dBenson et al.,,2003; .KJiochfar et alj 
l2007h . These semi-analytical models allow us to study morpho- 
logical evolution: how the morphology of a galaxy changes in 
time. To compare these models with the real Universe, we need 
to know the observed LF in detail. 

The influence of the global en vironment on the LF has been 
investigate d by H ovle et aH (l2005h . using the SDSS data, and by 
ICroton eral.l (2005), using the 2dFGRS data. These results show 
strong environmental trends: galaxies in higher density regions 
tend to be redder, of earlier type, have a lower star formation 
rate, and are more strongly clustered. Some of these trends can 
be explained with the well known m o rphology-density relation 
dEinasto et all Il974t iDressleH [l98^; 'Postman &GelleiJ 11984 
and luminosity-density relation CHamilton 19881) . It is less well 
known how far these trends extend when moving toward extreme 
environments, into deep voids or superclusters. 

Recent studies have shown that dust plays an impor- 
tant role in galaxy evolution and it may significantl y influ- 
ence the luminosities and colours of galaxies dPierini et alJ 
2004t iTuffs et al.1 12004 iDriver et alJl2007bt iRocha et alJl2008t 



Table 1. Galaxy parameters used. 



Tempel etalJbOlo C especially for late-type galaxies. Thus, in 



order to study intrinsic properties of galaxies, it is necessary 
to take dust exti nction into account. Using the SDSS data, 
IShao et al ] (l2007l) have studied the influence of dust on the LF. 
In general, dust is important for late-type spiral galaxies; nearly 
edge-on galaxies are most affected. 

In the present paper we use the SDSS data to study the LF in 
different global environments and for diff'erent types of galaxies, 
taking the effect of dust attenuation into account. The LF de- 
pendency on group properties will be analysed in a forthcoming 
paper 

Throughout this paper we assume a Friedmann-Robertson- 
Walker cosmological model with the total matter density Qn, - 
0.27, dark energy density Qa = 0.73, and the Hubble constant 
Ho = lOO/ikms 'Mpc '. Magnitudes are quoted in the AB sys- 
tem. 



2. Data and method of analysis 

2.1. Galaxy sample 

Our present analysis is based on the SDSS Data Release 7 
(lAbazaiian et al ]l2009h . We used the main galaxy sample, down- 
loaded from the Data Archive Server (DAS) of the SDSS. 
Additionally, some of the galaxy parameters that were not avail- 
able from the DAS, were downloaded from the Catalog Archive 
Server (CAS) of the SDSS using a Structured Query Language 
(SQL) search. We use only the contiguous imaging and spectro- 
scopic area of the Northern Galactic Cap (the Legacy Survey), 
which covers 7646 square degrees of the sky. The main galaxy 
sample of the SDSS includes galaxies brighter than the limit- 
ing r-band Petrosian magnitude 17.77. Our galaxy sample is de- 
scribed in detail in Tago et al. (2010); in total, this sample in- 
cludes 583362 galaxies. 

For the present work we use galaxies that lie in the distance 
inte rval 55-565 /i~'Mpc. W e use the co-moving distances (see, 
e.g. iMartmez & Saar"2002') as the linear dimensions, calculated 
using the cosmological parameters listed above. The lower dis- 
tance limit (55/i 'Mpc) is chosen to exclude galaxies of the 
Local Supercluster; the upper limit (565/i"'Mpc) has been set 



Name" 


SDSS Name* 


Description 


m/ 


- 


extinction-corrected 






cmodel magnitude 






absolute magnitude 


z 


z 


spectroscopic redshift 


d 




co-moving distance 


/deV 


fracDeVj-'' 


weight of the r''"* component 


?exp 


expAB_r 


exponential fit ajb 


9dcV 


deVABj 


de Vaucouleurs fit ajb 


'"exp 


expRadj 


exponential fit scale radius 



Notes. Parameter name, as used in the present paper. Parameter 
name as given in the SDSS CAS archive. x is either u, g, r, i or z filter 
in the SDSS; the magnitude is corrected for the Galactic extinction. 

Weight of the de Vaucouleurs component in the best-fit composite 
model. 

because the SDSS sample becomes very diluted at large dis- 
tances. Since the faint limiting magnitude has varied throughout 
the survey (between 17.5 and 17.77), we shall use the limiting 
magnitude 17.6 in the r-band; this leads to a uniform distribu- 
tion for the y/Vmax (see Sect. 12.2b . In addition to the lower limit- 
ing magnitude 17.6, we use the upper limiting r-band magnitude 
14.5, since the brighter end of the SDSS sample is incomplete 
due to saturation of bright galaxy images and due to blending 
with bright saturated stars. As suggested in the SDSS webpagqj, 
we use composite model magnitude^ 

All the galaxy parameters used are given in Table [1] together 
with their SDSS CAS archive names. The apparent magnitude m 
was transferred into the absolute magnitude M according to the 
usual formula 



= niA - 25 - 5 logio(t/z,) - K, 



(1) 



where the luminosity distance di - d(l + z); d is the co- 
moving distance in the units of /i 'Mpc and z is the observed 
redshift. The term K is the ^-He-correction. The apparent magni- 
tude m^i is the composite model ffux, calculated according to the 
SDSS webpage^, cor rected for Galactic extinction according to 
ISchlegel et an(ll998h . 

The ^-corrections for the SDSS galaxies were calculated us- 
ing the KCOR RECT alg orithm (version v4 1_4) d eveloped by 
iBlanton et all (l2003ah and lBlanton & RoweisI (l2007l). The evolu- 
tion cor rections e have been applied according to Blanton et al.l 
(l2003bl) . 



2.2. Selection effects and luminosity function estimation 

Several methods have been used in the past to estimate the lu- 
minosity functions: VmJv, C , various maximurn-hke lihood esti- 
mators (see, e.g. lWillnieill997l;IWall & Jenkinsl2003!) . The main 
reasons for so many methods were small sample volumes and, 
consequently, small numbers of galaxies. For such samples, de- 
viations from overall homogeneity inffuenced the results, espe- 
cially for the simplest, V^l^ method. 

Nowadays, both the sample volumes and sample sizes are 
about two orders of magnitude larger than in the recent past, and 
all methods work equally well. We chose the V^!,^ method as 
the physically most transparent and statistically straightforward. 
Since the luminosity function is practically a probability distri- 
bution, it should be estimated as a distribution, and this is what 



' http://www.sdss.Org/dr7/algorithms/photometry.html#which_mags 
" http://www.sdss.Org/dr7/algorithms/photometry.html#cmodel 
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the number density of galaxies of the luminosity L) was found 
as usual: 



250 300 350 
d (h'''Mpc) 

Fig. 1. The absolute magnitude of the SDSS DR7 main sample 
in the r-band versus distance. Black solid lines show the upper 
and lower limits of the flux-limited sample: 14.5 and 17.6 mag 
in the observer frame r-band, respectively, with the average K- 
corrections. 




150 200 



250 300 350 400 450 500 550 
d (h"'Mpc) 



Fig. 2. Upper panel: the weight factor Wd at different distances. 
Lower panel: the average global density in thin concentric shells 
and given in units of mean density (calculated using the weight 
factor Wd)- 



the Vj^^jj and the C methods do. The C method gives the in- 
tegral distribution, thus its errors are more difficult to estimate 
than for the V:^^,^ method. Estimating probability distributions 
is a well-developed topic in statistics, and we can choose all 
the modern tools - kernel densities, adaptive kernels, smoothed 
bootstrap forpointwise confidence intervals (see Appendix lAli. 

The only requirement for the V^^,^ method is that, on an av- 
erage, the spatial distribution of galaxies should be roughly uni- 
form. That can be easily tested; we did that. 

Below we give a short summary of the V^^^ method; a de- 
tailed description with our modifications is given in AppendixlAl 

In a flux-limited sample, the principal selection effect is the 
absence of galaxies fainter than the survey limiting magnitude. 
This eff'ect is well seen in Fig.[Tl showing the absolute luminosi- 
ties of galaxies plotted against their distance; at large distances 
only the brightest galaxies are seen. 

To take this effect into account in the determination of the 
LF of galaxies we used the standard V^),^ weighting procedure. 
The differential luminosity function «(L)dL (the expectation of 



«(L)dL = Yj 



Vmax(-L,) ' 



(2) 



where dL is the luminosity bin width, 1a(x) is an indicator func- 
tion, selecting galaxies that belong to a particular luminosity bin; 
VmaxiL) is the maximum volume where a galaxy of a luminosity 
L can be observed in the present survey; summation is made over 
all galaxies of the survey. This procedure is non-parametric, and 
gives both the form and the true normalisation of the LF. 

As said above, the spatial distribution of galaxies in the sam- 
ple should be uniform. This can be tested by examining the dis- 
tribution of y/Vmax, where V = V{d) is the part of the sam- 
ple volume closer th an the distance c/ to a galaxy (see, e.g. 
Wal l & Jenkinsl2003h - this distribution should be uniform. Due 
to the large size of our sample it is indeed practically uniform 
- the K-S test gives the probability (l-2)xl0"^ for the distri- 
butions to differ. The V/Vmax distribution is sensitive to sample 
incompleteness; choosing the limiting r-band magnitude 17.6 in- 
stead of the standard 17.77 was essential to make it practically 
uniform. 

For a further test, we generated several test samples of 
500000 galaxies within the same distance and magnitude limits 
as for the observational sample, using the Schechter luminosity 
function with typical parameters (a = -1.14, M* - -20.6), and 
our method recovered the initial parameters with deviations of 
0.03% {o-a = 0.0003, ctm. = 0.003). 

As the luminosity function is rapidly changing with lumi- 
nosity, especially at the bright end, varying bin widths should 
be used. This is most easily achieved by an adaptive kernel es- 
timation of the LF; the method is explained in the AppendixlAl 
Uncertainties of the LF (the confidence regions of the estimate) 
were estimated by the smoothed bootstrap method, also de- 
scribed in the AppendixlAl 

The maximum volume V,mx(L), where a galaxy of an ab- 
solute luminosity L could be observed is set by the minimum 
and maximum distance at which the corresponding apparent lu- 
minosity would fall within the luminosity limits. The distance 
limits are shown as solid lines in Fig.[T] These limits were calcu- 
lated, using Eq. ([T]i, where the apparent magnitude m was set to 
the survey limiting magnitudes: 14.5 and 17.6 in r-band. Since 
the A:-correction is colour- and distance-dependent, we used the 
colour-averaged values of the ^-correction for survey galaxies 
found as a function of distance. 



2.3. Determining environmental densities 

The environmental effects on the LF were studied at the level of 
the global density field, the so called supercluster-void network. 
We use the luminosity density for calculating the global density, 
assuming that it is proportional to the total matter density. 

In calculations of the luminosity density, selection effects 
typical to flux-limited samples had to be considered: in distant 
fields, only the brightest galaxies have been detected, while in 
nearby regions, the brightest galaxies are absent (see Fig. [TJ. 
To take this effect into account, we calculated the distance- 
dependent weight factor Wd 



Wd^ 



J^°°Ln(L)dL 
f^' Ln{L)dL 



(3) 
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where Li,2 = LqIO""**'*^©"'^' -^ are the luminosity limits of the 
observational window at a distance d, corresponding to the ab- 
solute magnitude limits of the wind ow Mi and M 2; we took 
Mq = 4.64 mag in the r-band (Bla nton & RoweisH2007h . The 
limits Ml and M2 have been chosen the same as for the LF (black 
lines in Fig.[TJ. In the latter equation, n(L) is taken to be the LF 
in the r-band for all galaxies (see Table|3]l, because we are inter- 
ested in the total luminosity density field. 

We assumed that every galaxy belongs to a group of galaxies, 
several of which may lie outside the observational window. The 
luminosity of the unobserved galaxies was taken into account 
by multiplying the observed galaxy luminosities by the weight 
Wd- Due to their peculiar velocities, the distances of galaxies are 
somewhat uncertain; if the galaxy belongs to a group, we use the 
group distance to determine the weight factor For the luminosity 
density field calculations we used the luminosities without the 
intrinsic absorption corrections. 

The luminosity function (calculated as described in the 
AppendixlAj was approximated by a double-power law: 



Table 2. Properties of different environmental regions. 



niL)dL oc (L/L'fil + (L/L*yy^-''^^^d(L/L*), 



(4) 



where a is the exponent at low luminosities (L/L*) 1, 6 
is the exponent at high luminosities (L/L*) » 1, y is a pa- 
rameter that determines the speed of the transition between 
the two power laws, and L* is the characteristic luminosity 
of the transition. A simil ar double-power law was also used 
by IVale & Ostrike^ (|2004|) to fit the mas s-luminosity relation 
in their subhalo model and bv ICoorav & Milosavlievic (2005) 
to fit the LF of central galaxies. Since several papers have 
sho wn that t he Schech ter function is not a best fit for the 
LF dBlanton e t al. 2005t iMercurio et al.ll2006l: lYang et al.ll2008t 
iTempel et al.l2009) . especially at the bright end, we shall use the 
double-power law for analytical approximation. Parameters for 
the LF approximations are given in Table |3] We give there also 
the fits for the Schechter function, to facilitate comparison with 
the results of other authors. 

We used the B3 spline kernel to calculate the density field. 
We chose a 8/i~'M pc scale smoothing for global densities 
dEinasto et al.ll2007 ). D etails of the density calculation are given 



Liivamagi et al.l (l2010i) . 



The upper panel of Fig. |2] shows the weight factor Wd as a 
function of distance; the lower panel of Fig. |2] shows the aver- 
age global density at diff'erent distances. When the weight factor 
is applied, the global density becomes roughly constant; the re- 
maining variations are due to the large-scale structure. 

We split our sample into four global density regions: voids, 
superclusters, and two intermediate regions. Both the void and 
the supercluster regions contain approximately 10% of all galax- 
ies. The rest of the galaxies are divided nearly equally between 
the two other regions. We designate these four regions as Dl, 
D2, D3, and D4, where Dl is the void region and D4 is the su- 
percluster region; D2 and D3 are the two intermediate density 
regions. The number of galaxies in each density region and the 
corresponding volume are given in Table |2] 



3. Dust attenuation in galaxies 

In this section we describe the necessary steps to correct the 
galaxy luminosities for internal attenuation. Since dust affects 
mostly late-type, spiral galaxies, we start this section with classi- 
fying galaxies into spirals and ellipticals. Since dust attenuation 
also depends on the galaxy inclination angle and galaxy colour 



Region 


Densities" 








V'' 


Dl 


0.0-0.8 


54942 


32963 


7836 


57.6 


D2 


0.8-2.0 


173320 


86239 


36718 


29.1 


D3 


2.0-5.0 


208210 


85749 


57965 


12.0 


D4 


5.0-00 


66884 


23123 


22026 


1.3 



Notes. Densities are given in units of the global mean density 
(0.01526 X 10'" /iLoMpc"'). <*> Percent of the total volume. 




Elliptical galaxies (visual classification') 
Spiral galaxies (visual classification) 



0.2 



0.4 
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0.8 



Fig. 3. Classification of galaxies into spirals and ellipticals. Red 
points (mostly in the upper right corner) are ellipticals and blue 
points (mostly below the dashed line) are spirals. Empty circles 
are elliptical galaxies classified by visual inspection; crosses are 
visually classified spiral galaxies. Grey points are galaxies where 
classification is unclear. The dashed line shows our spiral galaxy 
selection criterion (see text). 



(or galaxy type), we take the necessary steps to take that into ac- 
count. We will end this section by showing how our attenuation 
correction works. 



3.1. Galaxy classification 

Spiral galaxies have more dust than ellipticals and therefore 
the observed luminosities of spiral galaxies are more affected 
by dust. If we want to correct for dust attenuation, we have to 
know the morphology of a galaxy. Since dust affects mainly spi- 
ral galaxies, we will compose a sample of spiral and/or disc- 
dominated galaxies. We will use this sample to study the effect 
of dust attenuation. Additionally to spiral galaxies, we will com- 
pose also a sample of elliptical and/or bulge dominated SO galax- 
ies. We will use this elliptical galaxy sample only for compari- 
son, and not for a detailed study. 

As one source, we will us e the morpholog ical classifica- 
tion by the Galaxy Zoo project ("Linto tt et al.ll200 8). The Galaxy 
Zoo project has led to the morphological classification of nearly 
one million o bjects from the SDSS data by visual inspection. 
iBanerii et al.l (i2010) used the Galaxy Zoo data to develop a ma- 
chine learning algorithm (an artificial neural network) for galaxy 
classification; they have published also distribution histograms 
for different parameters of various types of objects: stars, spirals, 
ellipticals, and mergers. In this paper we use these histograms as 
the basis for selecting dominantly spiral or dominantly elliptical 
galaxies. Additionally, we will use the galaxy colour distribu- 
tions from'Lintott et alj ( 120081) . 

We will compare the distributions from the Galaxy Zoo 
project with our own visual classifications. We have classified 
a small sample (of nearly one thousand) galaxies in the Sloan 
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Fig. 4. Colour-magnitude diagram for spiral and elliptical galax- 
ies. Spirals are marked by blue dots (the lower cloud) and ellip- 
ticals by red dots (the upper cloud). 
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0.2 0.4 0.6 0.1 

Probability of being spiral/elliptical 

Fig. 5. Morpholo gical type prob ability distributions according 
to Huertas-Companv et alj (1201 ID of galaxies classified as spi- 
ral (blue soUd line) and elUpticals (red dashed line) in this paper. 
In our final sample, only galaxies with the given probability at 
least 0.5 were included. 



Great Wall region (lEinasto et al.l2010h . In this paper we use only 
galaxies with a clear classification to test the distributions taken 
from the Galaxy Zoo project. Additionally, we found another 
criterion (/dev) to classify galaxies into spirals and ellipticals. 

In the SDSS one of the main parameters that determines the 
type of a galaxy is the photometric parameter /jev (see Table[T]i - 
the point-spread function corrected indicator of galaxy morphol- 
ogy. The surface luminosity distribution of each galaxy in the 
SDSS has been fitted by the exponential and the de Vaucouleurs 
profiles. The best linear combination of these is used to represent 
the profile of the galaxy, and fd^y indicates the fr action of lumi- 
nosity contributed by the de Vaucouleurs profile. 'Bernardi et al.l 
(2005) used /dev > 0.8 to select early-type galaxies. Shao et al. 
(l2007h used /dev < 0.5 to select galaxies which are dominated 
by the exponential component (i.e. spiral galaxies). 

We combine the value /dev with the exponential profile axis 
ratio ^exp to do the primary classification. In Fig.[3]we show the 
/deV versus ^exp plot; crosses are spiral galaxies and circles are 
elliptical galaxies as classified bv lEinasto et alJ (l2010l) . As seen 
in this Figure, all galaxies, which have Q f^m < 0.4, a re spirals; 
the corresponding distribution in Banerii et alj (1201 Ol) confirms 
that. Ivan den Bosch & van d e Veni (|2009) have also shown that 
for elliptical galaxies, the axial ratio is mostly greater than 0.5. 
We also see that when moving toward lower values of /dev, the 
spiral dominated region becomes larger: the value of ^exp can 



be larger This is expected, since low values of /dev point to 
disc -dominated objects. In our classification we will take this be- 
haviour into account. The dashed line in Fig. [3] shows the limit 
for spiral galaxies. We will also leave unclassified those galaxies 
where /dev > 0.95 or /dev < 0.05, since classification at these 
extreme values may be rather uncertain. 

For all other galaxies that cannot be classified directly in 
the /dev versus ^exp plot, we use the galaxy colour as an ad- 
ditional parameter. Distributions of gal axy colours for diff erent 
type s of galax ies have been taken from lLintott et al.l (l2008 l) and 
Ban erii et al.l (2010); additionally, we have checked these distri- 
butions, using our own small sample of visually classified galax- 
ies. Rather tight (and thus secure) constraints have been chosen 
for the colour criteria, since each criterion is used to add galaxies 
to the sample. If a galaxy satisfied one criterion, then the follow- 
ing criteria were ignored. In the brackets we give the cumulative 
numbers of galaxies in the spiral/elliptical samples, after apply- 
ing this criterion. The criteria are: 

- all galaxies, where ^exp < 0.4 or ^dev < 0.4 (and 0.05 < 
/deV ^ 0.95), are classified as spirals dBanerii et alj|2010l) 
(60335/0), 

- galaxies, where ^exp < 0.9 - 0.8/dev (and 0.05 < /dev ^ 
0.95), are spirals (Fig.© (105819/0), 

- galaxies, where the rest-fra me colour M„-M, > 3.0, are el- 
lipticals (lLintottetal.ll2008[ and Fig.© (105819/40389), 

- galaxies, where the rest-frame colour M„-M,- < 1.9, are spi- 
rals (Lintott et al. 2008) (208483/40389), 

- gala xies, where the observer-frame nig-mr > 1 .2, are ellipti- 
cals (iBaneriietalJiJo iO) (208483/46237), 

- galaxies, where the observer-frame nig-nir < 0.6, are spirals 
(iBanerii et alJl201flh (216822/46237), 

- galaxies, where t/exp < 0.45 or o'dev <„0.45, are classified 
as spi rals dvan den Bosch & van de VenI 120091 ; iBanerii et al.l 
1201 Oh (248604/46237), 

- galaxies, where /dev > 0.7 and ^exp > 0.7, are ellipticals 
(Fig.© (248604/142023). 

The last criterion is included to add mostly ellipticals and/or 
bulge dominated SO galaxies into our elliptical sample. All 
the other galaxies, not classified according to the criteria 
given above, are used only to study the LF of the total 
galaxy sample. For the final refinement of the sample we 
exploited the automati c class ification recently published by 
Huertas- Companv et al.l ([2(inJ). In this paper, the authors as- 
sign for every galaxy a probability of being early- or late-type. 
Figure |5] shows the distribution of probabilities of being early- 
or late-type galaxies for our classified spirals and ellipticals. It is 
seen th at our classification agre es very well with the classifica- 
tion bv iHuertas-Companv et al.l ([2011). However, to further in- 
crease the reliability of our classification, we excluded from our 
sample all the galaxies, where this probability was less than 0.5, 
leaving us with the final sample of 228074 spirals and 124545 
ellipticals. 

Figure[3] shows our classification in the /dev versus ^exp plot: 
all galaxies, which are marked by blue points, are classified 
as spirals; all galaxies, which are marked with red points, are 
classified as ellipticals; grey points are non-classified galaxies. 
Figure |4] shows the classical colour-magnitude diagram for our 
galaxies; spirals are statistically bluer and fainter, and ellipticals 
are brighter and redder Since we use also other parameters for 
galaxy classification, some of the spirals are located in the re- 
gion, where mostly ellipticals reside: however, the number of 
these spirals (M„ - Mr > 2.3) is relatively small (20% of all 
spirals). Most of these red spirals are those where /dev is large: 
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Fig. 6. LF for different distance intervals. LF of spirals and all 
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scale, respectively. The filled areas show the 95% confidence re- 
gions of the LF. 
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Fig. 7. Apparent axial ratio versus apparent size for spiral galax- 
ies. Solid line is our fit for e dge-on spiral gala xies. Red points 
with errorbars are taken from I Shao et al.l (l2007h . 



All galaxies 
Spiral galaxies 
Elliptical galaxies 




0.9 



0.7 0.6 0.5 0.4 
Axial ratio Q^^p 



0.3 



0.1 



1 

^ 0.8 
o 

I 0.6 

CO 

a 0.4 

0.2 






All galaxies 




Spiral galaxies " 




Elliptical galaxies 









0.9 0.8 0.7 



0.6 0.5 

cos (/) 



0.4 0.3 0.2 



0.1 



Fig. 8. Upper panel: distribution of the apparent axis ratio for all, 
for spiral, and for elliptical galaxies. Lower panel: distribution of 
the restored intrinsic inclination angle for all, for spiral, and for 
elliptical galaxies. 



gle depends both on the visible and the intrinsic axial ratios, we 
can restore the intrinsic inclination angle only statistically. In or- 
der to consider dust attenuation, the inclination angle is needed 
only for spiral galaxies, but we will calculate it also for elliptical 
galaxies. 

Assuming that we know the axial ratios and the intrinsic in- 
clination angle of a galaxy, we can calcula te the apparen t axial 
ratio q = bla using the expression given bv lBinnevl(ll985h 



_ A + C - ^|{A-C)'^ + 
fl/ A + C+ ^J(A-C)^ + B^' 



(5) 



e.g. their luminosity profile is bulge dominated, but the visible 
axial ratio (g'exp) is small and therefore we have classified these 
galaxies as spirals. Practically all of the disc-dominated spirals 
are located below elliptical galaxies in this colour-magnit ude di- 
agram. Using the Galaxy Zoo data, Skib ba et"al] (l2009h found 
that red spirals account for nearly a quarter of all the spirals, 
which is in accordance with our findings. 

In our classification about half of the galaxies (45%) are spi- 
rals, about one quarter (25%) are ellipticals and for 30% of the 
galaxies, the classification is unclear 

To test the reliability of our classification, we calculated 
the luminosity functions for different distance intervals: 65-200, 
200-300, 300^00, and 400-565/;- 'Mpc. The LFs for diff^erent 
types of galaxies are shown in Fig.|6] This figure shows that the 
classification is statistically correct since the LF for all cases is 
nearly distance independent. 

3.2. Restoring the intrinsic inclination angle 

The intrinsic absorption is closely related to the morphology 
of galaxies. Additionally, attenuation depends on the inclination 
angle (Tuff^s et al. 2004; Driver et al. 2007b; Tempel et al. 2010). 
To take this eff'ect into account, we need to know the intrinsic 
inclination for every galaxy. Since the intrinsic inclination an- 



cos^ 6 1 f cos^0\ sin^ff 
^ - — (si" '*+^) + — ' (6) 

^ = (l-^j^cos0sin20, (7) 

In the last equation, 1 > ^ > ^, f is the ratio of the shortest semi- 
axis to the longest semi-axis and ^ is the ratio of the two longer 
semi-axes. is the inclination of a galaxy - the angle between 
the plane of the galaxy and the plane of the sky; is the angle 
between the longest semi-axis and the line of sight. 

The SDSS gives the apparent axial ratio q of each galaxy im- 
age. For spiral galaxies we use the r-band axial ratio <7exp, taken 
from the best fit of the image with an exponential profile con- 
volved with the point spread function (Stought on et al.ll2002l) . 

To calculate the inclination angle from Eq. (O, we use a 
statistical approach. We assume that the longest semi-axes are 
randomly oriented in space and the ratios ^ and ( are random, 
with different probability distributions for spirals and ellipticals. 
We use the Monte Carlo method to select random values for in- 
clination angle and axial ratios. 
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To find the intrinsic ratio of the shortest semi-axis to the 
longest semi-axis, we use galaxies with a small value of ^exp 
and assume that this value is the ratio ^. As demonstrated by 
IShao et aT] (l2007h . the apparent axial ratio depends on the appar- 
ent size and therefore, the value of ^ depends also on it. To find 
how the intrinsic thickness of an observed spiral galaxy depends 
on its apparent size, we plot the apparent axial ratio versus the 
apparent size (see Fig. |7]i and for every apparent size, we find 
the minimum thickness as shown with solid line in Fig. |7] For 
every spiral galaxy, we find the ratio ^ according to the scaling 
in Fig. |7] (solid line). We selected edge-on galaxies for deriv- 
ing the scaling, thus the ratio ^ should be s tatistically correct. 
In the same figure, results from IShao et aP (12007) are shown; 
our distribution remains somewhat lower than that of lShao et alj 
(12007). The difference arises from different methods for select- 
ing spiral galaxies. IShao et aTl ((2007) have used only /dev < 0.5 
to select spiral galaxies; this criterion alone is not reliable for 
smaller galaxies. There are many galaxies, where /dev > 0.5 but 
^exp < 0.5 - in our c lassification, these galaxies are spirals and 
are located below the' Shao et al.l ( 12007 *) region in Fig.Q 

For elliptical galaxies the ratio of the shortest semi-axis 
to the longest semi-axis is taken to be 0.7, with an 1-cr er- 
ror 0. 1 . This is in agreement with ^v an den Bosch & van de Veiil 
(l2009h . who derived the axial ratios for 13 elUptical galaxies. 
Additionally, the ratio ^ must be smaller than the apparent axial 
ratio by definition. 

Since the discs of spiral galaxies are not round an d ellip- 
tical galaxies are triaxial, the ratio ^ is less than one. iRvdenI 
( 200 4) showed that the intrinsic axial ra tio of discs is 0.85"^°o, 



which is slightly smaller than derived bv I Andersen & Bershadvl 
(l2002h (0.9!o "g). Both values are close to that derived for el- 
liptical galaxies bv lvan den Bosch & van de VenI ( |2009|) . In this 
paper we adopt the ^ distribution with the maximum at 0.9; to- 
ward larger values, the 2-cr deviation is 0.05, and toward smaller 
values this deviation is 0.1. Additionally we demand that the 
value of ^ has to be larger than the shortest-to-longest axis ratio. 
We will use the same ^ distribution both for spiral and elliptical 
galaxies. 

Figure [8] shows the distribution of q for all galaxies, for 
spiral, and for elliptical galaxies separately. If the ratio q were 
caused only by the inclination of a galaxy, then, assuming that 
the rotation axes were randomly oriented, the probability distri- 
bution of q should be flat. As we see, the distribution is not flat 
and decreases toward both smaller and larger axial ratios. 

The results of restoring the inclination angles are presented 
in the lower panel of Fig. [8] The distribution of inclination angles 
is not perfectly random, but is much closer to that, than the ap- 
parent axial ratio distribution, especially for elliptical and spirals 
separately. 



3.3. Accounting for dust attenuation in spiral galaxies 

In this section we describe the necessary steps to take dust at- 
tenuation into account. We will use the morphological classifi- 
cation, derived in Sect. 13.11 dust attenuation will be considered 
for spiral galaxies only. Late-type spiral galaxies have gener- 
ally more dust than early-type spirals, thus dust attenuation is 
higher there; we will take that additionally into account. It is 
known that blue galaxies have more dust and therefore the at- 
tenuation is larger; for redder galaxies dust attenuation is less 
impo rtant. The galaxy colo ur is also an indicator of the galaxy 
tvpe. lMunoz-Mateos et al.l (12009) have shown that dust attenua- 
tion for spiral galaxies has a large scatter and is nearly constant 
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Fig. 9. The dependence of attenuation on galaxy inclination for 
disc-dominated and bulge-dominated spiral galaxies. 
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Fig. 10. Luminosity functions for the edge-on and face-on sam- 
ples of spiral galaxies. Darker lines are attenuation-corrected 
LFs, dimmer lines are the LF for observed luminosities. Filled 
areas show the 95% confidence intervals. 
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Fig. 11. Attenuation-corrected LF for edge-on, face-on, and in- 
termediate inclination angle galaxies. LFs have shown for spi- 
rals, ellipticals, and all galaxies. The LFs of spirals and all galax- 
ies are shifted up by three and four units in logarithmic scale, for 
better separation of these three types. The filled areas show the 
95% confidence regions of the LF. 



for mid- and late-type spirals. For early-type spirals, the attenu- 
ation decreases. 

Additionally, dust attenuation also depends on the inclina- 
tion angle, as predi cted by Tuffs etjjL^004) and confirmed also 
observationally by iDriver et al.l (i2007bh. Recently , dust attenu- 
ation in M 31 was studied by Tempel et alJ (l2010l) : the derived 
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model allowed to investigate also the attenuation dependence on 
the inclination angle, as shown in Fig. |9] At lower inclination 
angles (nearly face-on oriented galaxies) the total extinction is 
low, since the line-of-sight optical depth of the dust disc is low. 
Attenuation increases when moving toward higher inclination 
angles and reaches its maximum for nearly edge-on galaxies. 
For exactly edge-on galaxies, dust attenuation decreases again, 
as the thickness of the dust disc is generally much lower than the 
thickness of optical components. Figure|9]shows that the dust at- 
tenuation - inclination angle relation depends also on the shape 
of the galaxy. For bulge dominated galaxies, the relation is flatter 
and attenuation is larger For disc dominated galaxies, attenua- 
tion is lower and increases more rapidly while moving toward 
edge-on orientation. 

The inclination-dependent attenuation curves, used in this 
paper are shown in Fig. |9] To classify galaxies into bulge- 
dominated or disc-dominated spirals, we use the parameter /dev: 
for lower values, the galaxy is disc-dominated and for larger val- 
ues, the galaxy is bulge-dominated. We use a linear interpolation 
in /deV when moving from disc-dominated galaxies to bulge- 
dominated galaxies. Since we do not know how the attenuation 
exactly depends on the inclination angle and galaxy colour, we 
have used only the general trends for correcting the LFs statis- 
tically. The calibration of the corrections has been carried out 
by comparing the LFs for galaxies at different inclination an- 
gles. For each galaxy subtype (colour) we modify the attenu- 
ation curves given in Fig. |9] to minimise differences between 
the shapes of the LFs for diff'erent inclination angles. For that, 
we multiply the attenuation curve by a factor x that depends on 
the galaxy rest-frame u-r colour For red spiral galaxies (with 
u-r > 2.2), X - 0.5 provided the best fit. This is expected, be- 
cause red galaxies tend to contain less dust. For galaxies with 
u-r < 1 .8 X = LO. For the intermediate galaxies, the factor x 
changes linearly with the u-r colour. 

Figure [10] shows the LFs for the observed and attenuation- 
corrected luminosities for the edge-on and the face-on samples 
of galaxies. For the edge-on sample cos 6 < 0.2, for the face- 
on sample cos 6 > 0.8. It is well seen that the observed LFs are 
quite different for the face-on and edge-on galaxies, while the 
attenuation-corrected LFs are quite similar. Figure[TT]shows the 
attenuation-corrected luminosity functions for the spirals and 
the ellipticals and for all galaxies together at different inclina- 
tion angles: edge-on, face-on, and intermediate inclination an- 
gles. The intermediate inclination angle sample is defined by 
0.5 < cos < 0.6. We see that the LF is nearly inclination inde- 
pendent. For elliptical galaxies, the LF is absolutely inclination 
independent. For the edge-on spirals and for all edge-on galaxies 
together, small differences are still noticeable. 

4. Results 

For all of the LFs presented in this section, the Schechter and 
double-power-law parameters are given in Table [3] 

4.1. Attenuation-corrected luminosity function 

Figure [12] shows the LFs for all galaxies, for spirals, and for 
ellipticals in the attenuation-corrected (lower panel) and uncor- 
rected (upper panel) case. At the bright end, most of the galax- 
ies are ellipticals, and at the faint end, most of the galaxies are 
spirals, as found in many previous studies. Correction for dust 
attenuation increases the brightness of spiral galaxies; the shift 
is especially noticeable at the bright-end of the LF. However, 
even with the attenuation correction applied, the brightest spiral 
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Fig. 12. LFs for all galaxies, spirals, and ellipticals. Upper 
panel: observed LFs. Lower panel: attenuation-corrected LFs. 
The filled areas show the 95% confidence regions of the LF. 
Light grey lines are analytical double-power-law functions. 



galaxies are still less luminous than the brightest elliptical galax- 
ies by about 0.5 mag. Table |3]gives the analytical fits to the LFs 
for the double-power law as well for the Schechter function. In 
Fig. [T2]the double-power-law fits are shown as solid lines for 
illustration. 

4.2. Luminosity functions in different environments 

Figure [T3] shows the LFs in different environments: Dl to D4, 
where Dl is the least dense (void) environment and D4 is the 
most dense (supercluster) environment. In each panel, lower and 
darker lines are the attenuation-corrected LFs and upper, dim- 
mer lines are the observed LFs. The most notable trend is that 
while moving from the lower global densities toward the higher 
ones, elliptical galaxies start to dominate the bright-end of the 
LF. In the least dense environments, ellipticals and spirals are 
equally abundant at the bright-end of the LF; in the densest en- 
vironments, the brightest galaxies are mostly ellipticals. At the 
faint-end of the LF, while moving from the low density regions 
to the high density regions, the difference between ellipticals 
and spirals decreases; in denser environments, the fraction of el- 
liptical galaxies increases. In general, differences between the 
attenuation-corrected and uncorrected cases are similar for all 
environments. 

A notable feature in the LF of ellipticals in the densest envi- 
ronments is a local minimum near Mr = -18 mag. A seemingly 
similar feature is present for spirals at M,- x -19 mag. However, 
these minima are of different origin: in the case of spirals, the 
small dip becomes visible because of an interplay with the bump 
at M,- a: -19.8 mag, mainly caused by a selection effect. The 
number of galaxies in the most dense environments is relatively 
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Fig. 13. LFs for different types of galaxies in different environments: upper left panel shows the least dense environment (Dl) and 
bottom right panel shows the most dense environment (D4). Green solid lines show the LFs for all galaxies; blue dashed lines show 
the LFs for spiral galaxies; red dotted lines show the LFs for elliptical galaxies. The upper LFs in each panel are the observed LFs 
that are shifted up two units in logarithmic scale and right by one mag. The lower LFs in each panel are the attenuation-corrected 
LFs. Filled areas show the 95% confidence regions. The panels Dl to D4 show different global environmental regions: Dl are void 
regions and D4 are supercluster regions. The LFs have been normalised to the volume of each sample. 



small in the SDSS sample volume and the presence or absence of 
a rich supercluster (the most dense environments) at a given dis- 
tance interval may leave a notable feature in the LFs because of 
the apparent luminosity limits of the survey. At M,. ^ -19 mag, 
the distance interval between 150 and 200h"'Mpc determines 
most of the LF (see Fig. [1). However, no rich superclusters are 
found in this region. 

In Fig. [14] the LFs for spiral galaxies, for elliptical galax- 
ies, and for all galaxies together are shown for different environ- 
ments. For elliptical galaxies, the bright-end of the LF moves to- 
ward higher luminosities when moving toward higher densities. 
This means that bright elliptical galaxies are residing mostly in 
high density environments, e.g., in the cores of galaxy clusters. 

Interestingly, the LF for spiral galaxies is almost indepen- 
dent of environment. The faint-end of the LF of the most dense 
environments is slightly different, but the number of galaxies in 
this region is also small and the dip of the LF may be caused 
purely by selection effects in the SDSS, as mentioned above. 
The bright-end of the LF of the least dense environments is also 
slightly different from that in other environments, because gen- 
erally, ver y bright galaxi es are absent from the low density envi- 
ronments dTempel et alJi200 9). 

Comparing the LFs of spiral galaxies and elliptical galaxies, 
the LFs for spirals slightly increase at the faint end in all envi- 
ronments, while the LFs of ellipticals have a maximum at about 
Mr ~ -20 and decrease toward lower luminosities. The LF for 
all galaxies is a combination of the spiral LF and the elliptical 



LF: it is determined by ellipticals at the bright end and by spirals 
at the faint end. 

Figure[T5]shows the LFs for red and blue ellipticals and spi- 
rals separately in different global environments. In general, the 
faint-end of the LF is mostly built up by bluer galaxies and the 
bright end includes mostly redder galaxies; this behaviour is the 
same for spirals and ellipticals. 

From Fig. [15] we also see that the increase of the number 
of bright ellipticals in dense environments is mostly caused by 
red ellipticals. As mentioned above, the LF of spirals is indepen- 
dent of the global environmental density (see Fig.[T4li. However, 
small changes with environment are seen at the faint-end of the 
LF. Figure[T5]shows that this change is much smaller for the red 
and blue spirals separately and is increased by the interplay of 
the differences in the LF shapes of these subpopulations. Once 
again, these differences concern only the densest environments 
and are thus a subject for selection effects. 

5. Discussion 

5.1. Interpretation of tiie results 

Our analysis of the LFs of galaxies of different morphology in 
different global environments shows that the global environment 
has played an important role in galaxy evolution. 

An interesting result of the present study is the finding that 
the LF of spiral galaxies is almost independent of the global en- 
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Fig. 14. LFs for different types of galaxies in different environ- 
ments. The upper LFs in each panel are the observed LFs that 
are shifted up two units in logarithmic scale and right by one 
mag. The lower LFs in each panel are the attenuation-corrected 
LFs. For elliptical galaxies (bottom panel), only the observed 
LFs are shown. Dotted lines stand for the least dense (void) envi- 
ronments; red solid lines stand for the most dense (supercluster) 
environments. The filled areas show the 95% confidence regions 
of the LFs. 



vironment, especially when looking at the red and blue spirals 
separately. However, one has to keep in mind that only the nor- 
malised distributions are similar; the actual number densities of 
spiral galaxies are different. On the contrary, for elliptical galax- 
ies, the global environment plays an important role: brighter el- 
liptical galaxies are located mostly in denser environments. In 
the case of ellipticals, the environment is more important for the 
red galaxies than for the blue ones. 

The results suggest that the evolution of spiral galaxies is 
slightly different for different types (colours) of spirals, but for 
a fixed-type spiral galaxy, the evolution is independent of the 
global environmental density. Thus the formation history of spi- 
ral galaxies in various global environments has to be similar. 
This result seems to be contradicting with the general ACDM 
cosmology: galaxy luminosity should be linked to the mass of 
the parent dark matter halo, the distributi on of which depends on 
the environment. Besides, according to iDelgado-Serrano et alj 
(1201 0). approximately half of the spirals were already in place 
6 Gyr ago and another half formed in mergers of irregular galax- 
ies. This suggests that the minor mergers and quiescent star for- 
mation are the dominant factors that determine the formation of 



spiral galaxies. A possible interpretation of our results may lie 
in the fragility of spiral galaxies: they form and survive only in 
specific conditions (e.g. the preservation of the gas, the absence 
of major mergers) which are typical of low density regions, but 
to some extent can be present also in high density regions. Some 
haloes may remain intact and host a spiral galaxy regardless of 
the global environment, while most of the potential hosts of spi- 
ral galaxies end up hosting ellipticals as demons trated by the 
semi-analytical models (Be nson & Devereuxll2010h . 

The derived LF of elliptical galaxies can be reconciled with 
the hierarchical galaxy formation through mergers. The denser 
the environment, the brighter galaxies there should reside be- 
cause of the increased merger rate. The difference between the 
LFs of elliptical galaxies in different environments is more no- 
table for red galaxies, in accordance with their supposed merger 
origin. This interpretation agrees well also with the picture of hi- 
erarchical formation of galaxies: for blue galaxies, the evolution 
is more quiescent and major mergers are not so important; for red 
ellipticals, merging is the dominant factor of galaxy evolution. 
Since blue ellipticals are most likely SO-s or late-type ellipticals, 
they have still some gas available for star formation and there- 
fore the evolution of blue ellipticals is closer to the evolution of 
spiral galaxies - the global environment is less important. 

5.2. Comparison with earlier work 

Influence of the global environment on the LF has been previ- 
ou sly studied by Hoyl e et al.l (l2005h using the SDSS data and 
by ICroton et al. (2005) using the 2dFGRS data. These works 
have shown that galaxies in the void environment are primarily 
of late-type. In our results, this effect becomes especially pro- 
nounced after applying the absorption corrections: Fig.fTSlshows 
that in void regions, spiral galaxies dominate over ellipticals at 
all luminosities. 

The LFs derived by iHovle et al.l (|2005|) for a considerably 
smaller sample of galaxies are quite similar in different envi- 
ronments, except for the hig hest density environ ment, where the 
faint-end slope is shallower. Croton et al.l ( l2005b have found that 
the faint-end of the LF depends weakly on environment. In gen- 
eral, our analysis confirms this result, except for the highest den- 
sity environment, where an excess of faint galaxies compared 
to other environments is found (noticeable for red spirals and 
for blue ellipticals in Fig. [TsT l. This excess has been detected 
also by Xia et aL (2006). In deep surveys of the Hubble Space 
Telescope, an exce ss of faint red ga l axies has been found in the 
field environments: ISalimbeni et al. (|2008|) hav e seen such trend 
in the GOODS dataset and lDrorv et alJ (l2009h in the COSMOS 
field. Thus the excess of faint red galaxies appears in all envi- 
ron ments, most s trongl y in dense cluster regions. 

iPhlepsetaPdloOTl) " have studied the global environment be- 
yond the redshift 2. Using three different fields with different 
global environments, they show that for blue galaxies, the envi- 
ronment plays a smaller role than for red galaxies. Their results 
are in agreement with our findings for the relatively nearby re- 
gion. 

Many previous works were concentrated only on the local 
environment. For example, Yang et al. (2009) studied the LF for 
the central and satellite galaxies in groups, using the SDSS data. 
They found that in general, red galaxies are the central galax- 
ies and blue galaxies are satellite galaxies; however, they found 
that for very low masses, the number of red central galaxies in- 
creases. They speculate that these galaxies are located close to 
large haloes so that their star formation is truncated by the large 
scale environment. Our results also show that the faint-end of 
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Fig. 15. Attenuation-corrected LFs for red {upper panels) and blue (lower panels); for spirals (left panels) and ellipticals (right pan- 
els), in different global environments. Dl is the least dense environment (voids); D4 is the most dense environment (superclusters). 
The filled areas show the 95% confidence regions of the LFs. 



the LP increases when moving toward very high density (see 
Fig.[T4ll. When splitting our galaxies into the red and blue sam- 
ples (see Fig. [ISl l. the increase of the LF at the faint end is no- 
ticeable for every subsample (see TableO, indicating a universal 
trend. In our case, the increase is also noticeable for blue galax- 
ies. Ho wever, a direct comparison with the results of lYang et al.l 
(l2009l) is difficult, because the environmental densities have been 
estimated differently. 

IZandivarez et (12006') shows that the local environment 
(galaxy group mass) is an important factor in galaxy evolution. 
They show that the faint-end slope is practically constant for the 
blue cloud galaxies, while for the red sequence galaxies, the faint 
end is steeper for more massive systems. Their results can be in- 
terpreted in terms of galaxy mergers as the main driving force 
behind galaxy evolution in groups. Using local environmental 
densities, we get similar results: for ellipticals (that are domi- 
nantly red), the faint-end slope is changing with density, and for 
spirals (that are dominantly blue), the faint-end slope is practi- 
cally constant. We shall discuss the group environment in more 
details in a forthc oming study, using the group catalogue from 
iTago et al.l(l2010h . 

In general, the characteristic magnitudes for ellipticals are 
brighter th an for spirals a nd the faint-end slopes are steeper 
for spirals. iDevereux et al.l (12009,) used the /T-band luminosity 
to derive the LF for different Hubble types; classification was 
performed visually. The average shape of the LFs of ellipticals 
and spirals is generally in agreement with our results. In ad- 
dition, they found that the faint-end slope is steeper for late- 
type spirals than for early-type (SO) spirals and that the char- 
acteristic luminosity is larger for early-type galaxies. In Fig. [16] 
we use the value /dev to separate spiral galaxies into early 
an d late type. This figure shows a similar trend as pointed out 
by IDevereux et al.l ( l2009l) . but in our case, the differences are 
smaller. The qualitative results presented in Fig. [T6]remain the 




-17 -18 -19 -20 -21 -22 -23 

Mf- 5log/i (mag) 



Fig. 16. The attenuation-corrected LF for spiral galaxies for var- 
ious Hubble types. The Hubble type indicator is the value of /dev^ 
low value - late-type; high value - early-type. Errors are omitted 
for clarity. The filled areas show the 95% confidence regions of 
the LF 

same, when using the classification by iHuertas-Companv et al.l 
(I20T1I) . 

6. Conclusions 

We have used the SDSS data to derive the LF of spiral galaxies, 
elliptical galaxies, and all galaxies together in various environ- 
ments. We have taken special care to correct the galaxy lumi- 
nosities for the intrinsic attenuation. The principal results of our 
study are the following: 

- The LF of elliptical galaxies depends strongly on the envi- 
ronment; this suggests that global environmental density is 
an important driving force (via merging history) of ellipti- 
cal galaxy formation. Density environment is more impor- 
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tant for red elliptical galaxies than for blue elliptical galax- 
ies. 

- The evolution of spiral galaxies (the LF of spiral galaxies) is 
almost independent of environment, especially for blue and 
red spirals separately, showing that spiral galaxy formation 
has to be similar regardless of the surrounding global density. 

- The highest global density regions (superclusters) are sig- 
nificantly dif ferent from other regions, as stated also by 
iTempel et al.l 12009). Notably more faint galaxies are found 
in high density regions than in other environments. 

- The brightest galaxies are absent from the void regions. After 
correcting for the intrinsic absorption, spiral galaxies domi- 
nate the LF of void regions at every luminosity. 

- The faint-end of the LF is determined by spiral galaxies and 
the bright end by elliptical galaxies. The faint end includes 
mostly blue galaxies and the bright end mostly red galaxies. 

- Detailed studies of LFs require galaxy luminosities to be cor- 
rected for the intrinsic absorption by dust. Dust absorption 
affects mostly the bright-end of the LF. For the full LF, in- 
cluding all galaxies, the characteristic luminosity increases 
after attenuation correction. The faint-end slope of the LF is 
practically independent on dust attenuation. 

A comparison of these results with predictions of numerical 
simulations and/or semianalytical models would provide strin- 
gent constraints on the driving factors of the formation and evo- 
lution of galaxies in dark matter haloes. 
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Table 3. The parameters for the Schechter and double-power law LFs. 



Schechter 



Sample 



Double-power law 
T 6 



M' 



All galaxies (observed) 
Spiral galaxies (observed) 
Elliptical galaxies (observed) 


-1.12 ±0.01 
-1.18 ±0.01 
-0.19 ±0.01 


-20.71 H 
-20.17 H 
-20.59 H 


tO.Ol 
tO.Ol 
tO.Ol 


-1.30 H 
-1.42 H 
-0.26 H 


tO.Ol 
tO.Ol 
t0.02 


1.73 ±0.05 

1.74 ±0.05 
1.24 ±0.03 


-7.78: 
-9.52 : 

-14.64 


tO.36 
tO.52 
± 1.27 


-21.98 ±0.06 
-21.69 ±0.07 
-22.96 ±0.14 


All galaxies (Dl, observed) 


-0 


98 H 


to 


02 


-20 


04 H 


tO 


01 


-1 


36 H 


tO 


01 


2 


00 = 


tO 


05 


-12.70 


± 1.13 


-21 


81 


± 





08 


All galaxies (D2, observed) 


-0 


98 H 


to 


02 


-20 


42 H 


tO 


01 


-1 


32 H 


tO 


01 


1 


79 = 


tO 


03 
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± 1.21 


-22 


36 


± 





08 


All galaxies (D3, observed) 
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01 H 


tO 


01 
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tO 


01 


-1 


24 H 


tO 


01 


1 
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tO 


06 
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± 1.11 


-22 
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12 


All galaxies (D4, observed) 
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tO 


02 
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tO 
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-1 


18 H 


tO 


05 


1 
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tO 
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58 
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50 


Spiral galaxies (Dl, observed) 
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02 
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All galaxies (dust-corrected) 
Spiral galaxies (dust-corrected) 
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-1.15 


±0.01 
±0.01 


-20.66 ± 
-20.32 ± 


0.01 
0.01 


-1.29 H 

-1.41 H 


tO.Ol 
tO.Ol 
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1.79 
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-21.88 ±0.06 


All galaxies (Dl, dust-corrected) 
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-1 
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0.26 
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t0.05 
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-20.02 = 
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± 


0.27 
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±0.06 


-14.40 = 


t2.00 


-22.41 


± 
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±0.15 
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± 


0.16 


Blue spirals (Dl, dust-corrected) 
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0.01 


-1 
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t0.03 
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± 
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± 
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-20.19 ± 


0.01 


-0 


52 H 


t0.04 


1.34 


±0.08 


-20.00 = 


t2.00 


-22.82 


± 


0.49 


Blue ellipticals (D3, dust-corrected) 


-0.48 


±0.01 


-20.48 ± 


0.01 


-0 


63 H 


t0.03 


1.27 


±0.07 


-20.00 = 


t2.00 


-23.15 


± 


0.57 


Blue ellipticals (D4, dust-corrected) 


-0.88 


±0.02 


-20.86 ± 


0.02 


-0 


61 H 


t0.12 


0.85 


±0.16 


-20.00 = 


t2.00 


-24.16 


± 


0.50 



Notes. M* is in units of mag-5 \og{h). 
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where dL is the luminosity bin width, lAix) is the indicator func- 
tion that selects the galaxies belonging to a particular luminosity 
bin, VmaxiL) is the maximum volume where a galaxy of a lumi- 
nosity L can be observed in the present survey, and the sum is 
over all galaxies of the survey. 

This approach gives us the binned density histogram that de- 
pends both on the bin widths and the locations o f the bin edges; 
a better way is to use kernel smoothing (see, e.g. lWand & JonesI 



Appendix A: Confidence intervals for luminosity 
functions 

The standard V^l^ weighting procedure finds the differential lu- 
minosity function n(L)dL (the expectation of the number density 
of galaxies of the luminosity L) as foUows 
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Il995h . where the density is represented by a sum of kernels cen- 
tred at the data points: 

The kernels K{x) are distributions (_K{x) > 0, J K(x)dx = 1) of 
zero mean and of a typical width h. The width h is an analogue 
of the bin width, but there are no bin edges to worry about. 

As the luminosity function is rapidly changing with luminos- 
ity, especially at the bright-end of the LF, the bin widths should 
vary. This is most easy to implement by adaptive kernel estima- 
tion of the LF - instead oflA.2lwe write 



^ Vmax(i;) hi \ hi j 

where the kernel widths depend on the data, hi - h(Li). 

The choice of the kernel widths is a matter of ongoing study , 
but recommendations are available (see, e.g. I Silvermanll 1 997h . 
The kernel widths are known to depend on the density f(x) itself, 
with h ~ /(x)"'^^ for densities similar to normal distribution. 
This choice requests a pilot estimate for the density that can be 
found using a constant width kernel. 

We used the magnitude scale for our luminosity function (all 
kernel widths are in magnitudes), and the Bj box spline kernel: 

Bi(x) = (|x-2|^ -4|x- + 6|x|^- 
-4|x+ H- |jt;H-2|^)/12. (A.4) 

This kernel is well suited for estimating densities - it is com- 
pact, differing from zero only in the interval x € [-2, 2], and it 
conserves mass: 2,- B^ix - /) = 1 for any x. 

For the pilot estimate, we used a wide kernel with the scale 
h = 0.5 mag. For the adaptive kernel widths, we adopted h - 
0.05 mag (the typical SDSS rms magnitude error) as the minimal 
width (for the maximum pilot density) and rescaled it by the 
h ~ /piiot(jc) '^^ law. The luminosity function drops sharply at 
the bright end, leading to very wide kernels; we restricted the 
kernel width by 2/i - 0.5 mag. 

If we choose the kernel width this way, we minimise the 
mean integrated standard error (MISE) of the density. While that 
is certainly a useful quantity, we are also interested in the "error 
bars", pointwise confidence intervals for the density. These can 



be obtained by smoothed bootstrap (Silverman & Younji 19871: 
iDavison & Hinklevlll997tlFioridl2004r Here the data points for 
the bootstrap realisations are chosen, as usual, randomly from 
the observed data with replacement, but they have an additional 
smoothing component: 

L* = Lj + hsj, (A.5) 

where e is a random variable of the density K{x). 

We generated 10000 bootstrap realisations, using the adap- 
tive kernel widths as for the true luminosity function estimate. 
We show the centred 95% confidence regions in our figures. 
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